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ABSTRACT 

The NRAO VLA Sky Survey (NVSS) is the only dataset that allows an accurate determination of 
the auto-correlation function (ACF) on angular scales of several degrees for active galactic nuclei at 
z ~ 1. Surprisingly, the ACF is found to be positive on large scales while, in the framework of the 
standard hierarchical clustering scenario with Gaussian primordial perturbations it should be negative 
for a redshift-independent effective halo mass of order of that found for optically selected quasars. We 
show that a small primordial non-Gaussianity can add sufficient power on very large scales to account 
for the observed NVSS ACF. The best-fit value of the parameter /nl, quantifying the amplitude of 
primordial non-Gaussianity of local type, is /nl = 62 ± 27 (1 ct error bar) and 25 < /nl < 117 (2 a 
confidence level), corresponding to a detection of non-Gaussianity significant at the ~ 3(j confidence 
level. The minimal halo mass of NVSS sources is found to be M„ii„ = 10i2.47±o.26/j-i^^ (j^^) 
strikingly close to that of optically selected quasars. We discuss caveats and possible physical and 
systematic effects that can impact on the results. 

Subject headings: cosmological parameters - cosmology: theory - galaxies: halos - large-scale structure 
of universe 



1. INTRODUCTION 

The investigation of primordial non-Gaussianity offers 
a powerful way of testing the generation mechanism of 
cosmological perturbations in the early universe. Al- 
though the standard single-field, slow-roll, canonical ki- 
netic energy and adiabatic vacuum state inflation gener- 
ate very small non-Gaussianity, any inflationary model 
that deviates from t his may entail a large r level of it 
, (|Bartolo et all 120041: iKomatsu et al. 1 120091 and refer- 
ences therein). 

Deviations from Gaussian initial conditions are com- 
monly taken to be of the so-called local type and param- 
eterized by the dimensionless parameter /nl: 

; $-'/' + /nl(0'-(0')) , (1) 

where $ denotes Bardeen's gauge-invariant potential and 
\ </) is a Gaussian random fleld. In this Letter, we use the 
cosmic microwave background (CMB) convention for the 

quoted /nl values. 

; A method (Da lai et al. I l2008t iMatarrese fc Verde I 
, l2008f l for constraining non-Gaussianity from large scale 
structure (LSS) surveys exploits the fact that the clus- 
tering of extrema (i.e., dark matter halos where galax- 
ies form) on large scales increases (decreases) for pos- 
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itive (negative) /nl- In particular, a non-Gaussianity 
described by Equation (jlj introduces a scale-dependent 
boost of the halo power spectrum proportional to 
on large scales (fc < 0.03/i/Mpc), which evolves roughly 
as (1 -I- z). 

Extragalactic radio sources are uniquely well suited to 
probe clustering on the largest scales: (1) radio surveys 
are unaffected by dust extinction which may introduce 
spurious features reflecting the inhomogeneous extinc- 
tion due to Galactic dust; (2) due to their strong cos- 
mological evolution, radio sources are very rare locally, 
so that radio samples are free from the profusion of lo- 
cal objects that dominate optically selected galaxy sam- 
ples and tend to swamp very large-scale structures at 
cosmological distances; (3) thanks to the strong cosmo- 
logical evolution, even shallow radio surveys reach out 
to subst antial redshifts. The NRAO VLA Sky Survey 
(NVSS; iCondon et aTl 119981) offers the most extensive 
sky coverage (82% of the sky to a completeness limit 
of about 3 mJy at 1.4 GHz) with sufficient statistics to 
allow an accurate determination of the auto-correlation 
function (ACF), w{9\ on scales of up to several degrees 
(|Blake fc Wallll200a ICVerzier et al.l[200l . 

When a realistic redshift distribution of the NVSS 
sources is adopted, the interpretation of the measured 
w{9) in the framework of the standard hierarchical clus- 
tering scenario with Gaussian primordial perturbations 
requires an evolution of the bias factor radically differ- 
ent from that of optically sel ected QSOs (jNegrello et al. I 
120061 : iMassardi et al. II2O1O0 . in stark contrast with the 
similar evolution of the luminosity function. In fact, the 
observed w{9) is positive up to large (~ 10°) angular 
scales, which, for the median source redshift (z^ — 1), 
correspond to linear scales wher e the correlation function 
should be negative (see also iHernandez-Monteagudo I 
(2009)). Here, we explore whether the (scale-dependent) 
large-scale non-Gaussian halo bias could reproduce the 
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observed shape of the NVSS sources ACF, preserving 
the kinship with opticahy selected active galactic nuclei 
(AGNs). 

2. NVSS AUTO-CORRELATION FUNCTION 

We include in our analysis only NVSS sources brighter 
than 10 mJy, excluding the strip at |6| < 5°, where 
the catalog may be substantially affected by Galac- 
tic emissions. Th is ensures a uniform sky coverage 
(jBlake fc Wallll2002[ ): the effect of possible residual large- 
scale gradients producing and offset of the ACF is dis- 
cussed in §131 The NVSS source surface density at this 
threshold is 16.9 deg~^. The rcdshift distribution has re- 
cently been determined by [Brookes ct al. (2008). Their 
sample, complete to a flux density of 7.2 mJy, comprises 
110 sources with S'i.4ghz > lOmJy, of which 78 (71%) 
have spectroscopic redshifts, 23 have redshift estimates 
via the K-z relation for radio sources, and 9 were not 
detected in the K band and therefore have only a lower 
lir nit in z. Here, we have adopted the description given 
bv lde Zotti eTahl (pOlO l: 

dN/dz = 1.29 + 32.37z-32.89z^ + n.l3z^-1.25z'^ . (2) 

The NVSS m aps are pixelized usi ng the HEALPix soft- 
ware package (jCorski et al. l[2005l l with A^side = 64, cor- 
responding to A^pix = 49, 152 pixels with dimensions 
0.92° X 0.92°. 
The ACF estimator 'w{6) reads 



(n,; - fin){nj - fju) 



(3) 



where fi and are the coverage fraction and the num- 
ber of radio sources in each pixel, respectively; h is the 
exp ectation value for the number of objects in the pixel 
(see lXia et aFI [20091) . The sum runs over all the pixels 
with a given angular separation 9. The equal weight- 
ing used here is nearly optimal because of the uniform 
NVSS sky coverage and because on large scales the noise 
is dominated by sample variance. For each angular bin 
centered around 6*, Ng is the number of pixel pairs sepa- 
rated by an angle within the bin, weighted by the cover- 
age fractions. We used A^b = 9 angular bins in the range 
1° < < 8° with a linear binning plus another estimate 
at 9 = 40', since below 40' the correlation function is af- 
fected by multiple source components and above 8° the 
signal may be affected or even domi nated by spurious 
density gradients (.Blake fc Wallll2002D . 

We estimated the covariance matrix o f the data points 
using the jackknife re-sampling method fScra nton et al. 1 
[2002). We divide the data into M — 30 patches, then cre- 
ate M subsamples by neglecting each patch in turn, and 
in each subsample we measure the ACF. From the M 
estimates of the ACF functions, we estimate the diago- 
nal (variance) and off-diagonal (covariance) elements of 
the covariance matrix for w{9). In Figure [1] we plot the 
observed NVSS ACF which is consistent with previous 
estim ates using different approaches (e.g. I Blake fc Walll 
I2OOI . For comparison, we also show the best fit theoret- 
ical ACF curve for the Gaussian case (/nl = 0) assum- 
ing the redshift- independent minimal halo mass Mmin. 
Clearly, the curve does not match the observed ACF data 
on large scales. 
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Fig. 1. — Observed ACF of NVSS catalog. Values are jackknife 
estimated. The black solid line is the best fit model of our non- 
Gaussian calculations, while the red dashed line refers to the Gaus- 
sian case. The vertical arrow marks the angular scale above which 
the theoretical Gaussian ACF becomes negative. (Here, negative 
values are not visible, due to their very small amplitudes. See 
Figure[2]for details.) 



Note that the integral over the full survey solid angle 
covered by the observationally determined ACF vanishes 
by construction. The estimated values go to zero for 
9 ~ 30°, and the non-Gaussian model ACF shown in 
Figure [1] becomes negative approximately at the same 
9. However, no special meaning should be attached to 
this coincidence since, as noted above, the observational 
estimate of w{9) are unreliable for 9 > 8°. 

3. METHOD 
3.1. Modeling the effects of non-Gaussianity 

The effects of non-Gaussianity on the source clustering 
properties arise because a non-zero /nl affects the halo 
mass function and enhances the halo clustering on large 
scales. The second effect is the dominant one. 

In the presence of non-Gaussianity, the mass function 
riNcl-^: /nl) can be written in terms of the Gaus- 
sian one nQ™(M, z), for which a good fit to the results 
of simulations is provide d by the Sheth-Tormen formula 
(jSheth fc Tormeii1ll999D. multipl ied by a non-Gaussian 
corre ction factor (jMatarrese et al. .2000.:.LoVerde et al. I 
[200I) : 



i?NG(M, 2:,/nl; 



'M 



6(5cc(-z) 



S: 



3,M 



'M 



dS: 



3,M 



,(4) 



where the normalized skewness of the density field 
Ss^M oc /nli and (Tm denotes the rms of the dark mat- 
ter density field linearly extrapolated to z = and 
smoothed on the scale R corresponding to a Lagrangian 
radius of a halo of mass M. Here, 6cc{z) denotes the 
critical density for ellipsoidal collapse, which for high 
peaks is Sedz) ~ Sc{z)yi^ {q = 0.75) and has been cali- 
brated on N-body simulations (jGrossi et al. I [2009.) and 
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Fig. 2. — Effect of different Mmin masses on the ACF: the zero- 
crossing angular scale of the ACF decreases with increasing Mmin 
for /nl = 0. 



Sc{z) — Ac{z)D{0) / D{z) where D{z) denotes the hnear 
growth factor; Ac{z) ^ 1.68 and evolves very weakly 
with redshift. 

More importantly, the large-scale halo bia s is also mod- 
ified by the presence of non-Gaussianity (Dalai et al. I 
[2008; Matarrese k Verde 2008; Grossi et al. 200^7 

b-Nciz) - boiz) ~ 2(60(2) " l)f-NhScciz)aM{k) , (5) 

where the factor auik) encloses the scale and halo mass 
dependence. In practice, we find that, on large scales, 
ayiik) cx and is independent of the halo mass. 

We start by assuming that the large-scale, linear h alo 
bias for the Gaussian case is (ISheth fc Tormen lll999l) 



6g = 1 



1 



D{z,) 
2p 



qScizf) 



1 



dciz{)D{Zo) 



'M 



1 



(5c (zf). 
qS^{zi) 



1 P 



(6) 



where Z{ is the halo formation redshift, and Zo is the 
halo observation redshift. As we are interested in mas- 
sive halos, we expect that Zf ~ Zq. Here, q = 0.75 and 
p = 0.3 account for non-spherical collapse and are a fit 
to num erical simulations (see also [Mo fc White 19961: 
IMo. Jing fc White iri997tlScoccimarro et al. 1120011) . We 
will later relax this assumption. 
Finally, the weighted effective halo bias is given by 



,z,k,ff^i^) = 



(7) 



Two things should be clear from Equation ([S]): (1) 
there is a degeneracy between 6g and /nl (the same 
amount of non-Gaussian bias can be given by differ- 
ent pairs of Bq, /nl values; strictly speaking, is not 
a free parameter here, and the degeneracy is between 
-^^min: which is a free parameter, and /nl! however be is 
strongly dependent on Mmin); (2) the 1/k^ scale depen- 
dence means that large-scales are mostly affected by /nl 
and small scales are primarily affected by 6g ■ A positive 
/nl enhances the amplitude of auto-correlation power 
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Fig. 3. — Effects of non-Gaussianity on the auto-correlation 
power spectra (left panel) and on the ACF (right panel) for three 
different Gaussian bias models. 



spectra especially at large angular scales {£ < 200, 9 > 
4°). This is the effect we shall use to constrain /nl and 
its impact on the ACF is clearly visible in Figure [T] In 
fact, in the Gaussian case, for the adopted redshift dis- 
tribution and a redshift-independent Mmin, the ACF is 
expected to become negative for 6* > 4°. This is also 
shown in Figure [5] where the ACF is plotted for several 
values of Mmin. 

In general, the expected degeneracy between /nl and 
6g may be lifted in two ways: (1) on small scales the 
effect of /nl is completely negligible but not that of 6g, 
(2) the redshift dependences of the two contributions are 
different. However, the angular correlation function en- 
compasses the signal from different redshifts and differ- 
ent physical scales, complicating the separation of the 
/nl and be signals. 

To explore the effect of relaxing the assumption Zf ~ 
Zo, we have a lso considered a mode l for the Gaussian 
bias given by ([Matarrese et al. II1997I : iMoscardini et al. 1 
119981 and references therein) 

^g(2)=6i + -^ (0<7<2), (8) 

with bi and 62 being free parameters. Indeed, an "object- 
conserving" bias model corresponds to 7 ~ 1, while 
the bias of high-density peaks for objects that have just 
formed yields 7^2. In Figure [3l we show the effect 
of a Gaussian bias model on the auto-correlation power 
spectra and ACF when varying the power-law index 7; 
the larger the value of 7 the larger is the large-scale non- 
Gaussian boost. 

iLin I ([2001[ ) and lYoo. Fitzpatrick fc Zaldarriaga I 
([20091 ) discussed the gauge dependence of matter 
power-spectrum on very large scales (fc < 0.003 h/Mpc). 
We find that this gauge-dependent effect on the matter 
power spectrum can be mimicked by that of a non- 
Gaussian halo bias model with /nl ~ 5. Here, we 
calculate the matter power spectrum in the conformal 
Newtonian gauge. 



4 



Xia et al. 



3.2. Implementation and data sets 

The theoretical prediction for the ACF depends on: the 
cosmological parameters, the minimal halo mass Afmin 
and /nl- For the generalized bias model of Equation 
([5]), we also add the hi and 62 and 7 bias parameters. 

Rather than fixing the cosmological parameters to 
the best-fit values derived from Wilkinson Microwave 
Anisotropy Probe seven year (WMAP7), we will report 
the results after having marginalized over them with a 
prior given by a compilation of recent data sets. 

We perform a g lobal fitting using the CosmoMC package 
(jLewis fc Bridle I I2002I) . a Markov Chain Monte Carlo 
code, which has been modified to calculate the theo- 
retical ACF. We assume purely adiabatic initial con- 
ditions and a flat universe, with no tensor contribu- 
tion. We vary the following cosmological parameters 
(r^b/i^, ^ch"^, T, Qs, ns, As, /nl, A^min), where r^b/i^ and 
Qch'^ are the baryon and cold dark matter densities, t is 
the optical depth to reionization, Qs is the ratio (multi- 
plied by 100) of the sound horizon at decoupling to the 
angular diameter distance to the last scattering surface, 
Hs is the primordial spectral index, and is the primor- 
dial amplitude. We do not consider massive neutrinos 
and dynamical dark energy and for the pivot scale we 
set fcso = 0.05 Mpc~^ 

We also use (1) CMB temperature and polariza- 
tion angular power sp ectra as measured by WMAP7 
([Komatsu et al. I I2010D . (2) baryonic acoustic oscilla- 
tions in the galaxy power spectra as measured by the 
SDSS7 and the Two-degree Field Galaxy Redshift Sur- 
vey (2dFGRS; Percival et al. 2010) (3) SNIa distance 
moduli of Union compilati on from the Supernova Cos- 
mology Project (|Kowalski et al. 1120081 ). We add a prior 
on the Hubble constant, Hq = 74.2 ± 3.6 km/s/Mpc 

less et aL~l l2009t ) . Finally, we set the minimal halo 
mass Minin > 10^^h~^MQ consistent with observations 
showing that radio AGNs are hosted by halos more 
mass ive than those hosting optical QSOs (Hickox et al.l 
[2OOI). 

4. RESULTS 

We start by considering the case where the Gaussian 
bias is given by Equation (|6]). In Figure |4l we show 
the one-dimensional posterior probability distributions 
for Mmin and /nl after marginalizing over the other pa- 
rameters. The right panel shows the degeneracy between 
Mmin and /nl- Note that the constraints on /nl only 
come from the ACF data; external data sets are only 
used to set the underlying cosmology. 

We find that the current ACF implies /nl > at 
the ~ 3 cr confidence level. The reason for that can 
be clearly seen in Figure [2] where the Gaussian model 
with the redshift-independent My^in yields an ACF drop- 
ping to zero for ~ 4° and becoming negative on larger 
scales where the observed ACF is still positive. Non- 
Gaussianity of the local type adds power on large angular 
scales yielding a good fit to the observed data points. 

The marginalized constraints on the non-Gaussianity 
parameter /nl, 

/nl = 62 ± 27 (1 a CL) , (9) 
(6)25 < /nl < 117(142) [95% (99.7%) CL], (10) 




Fig. 4. — Marginalized one-dimensional and two-dimensional dis- 
tributions (1, 2(7 contours) of the minimal halo mass A/niin and of 
the non-Gaussian parameter /nl- 



are compatible with other previous esti mates 

(lYadav fc Wandelt I [2008t ISlosar et al. I [20081: 

iCurto et al. "2009';'Smidt et al. "2009'; 'Jimenez & \^rdej 
2009; Smith et al. 2009; Rudjord et al. 2003) and in 
very good agreement w ith the recent WMAP7 estimate 
(jKomatsu et al. Il2010t ). 

The minimal effective halo mass, Mmin — 
lQi2A7±o.26f^-ij^^ (Icr), turns out to be remark- 
ably close to that found for opti cally selected QS Os: 
Mqso = (3.0 ± 1.6) X lOi2/i-iM0 (jCroom et al.ll2005[ ). 

To explore whether a more general bias model (with 
7 allowed to conservatively vary even in an un-physical 
range 7 < 1) than that of Equation ^ may reconcile a 
Gaussian model with the data, we repeated the analysis 
using the bias of Equation ([8]). We keep bi and 62 fixed 
to 1.1 and 0.6, respectively, and we vary < 7 < 2 
(the values for bi and 62 have been chosen in order to 
provide a good fit to the data). In this case, the recovered 
central value of /nl becomes a little higher and the 1 a 
error bar increases by about a factor of 2 (a smaller /nl 
accommodates larger 7 and larger bias). 

We also perform a cross-check by using the published 
version of the ACF of Blake & Wall (2002) (in this 
case zero covariance between data points is assumed, 
since no covariance matrix has been computed) and find 
/nl = 58 ± 12 (Icr CL) (we rely on 12 data points in 
the range 0.5623° < 6* < 7.70795°). If we instead set 
the off-diagonal terms of our ACF estimate to zero, we 
get /nl = 70 ± 15 (1 cr CL). From these results we can 
conclude that covariance between data points increases 
the error bar on /n l by almost a factor of 2 and the 
IBlake fc Wall ((200^ and our ACF measurements are in 
very good agreement with each other in terms of de- 
rived /nl values. Another instructive cross-check is to 
see to what extent our conclusions are affected by apply- 
ing an overall subtraction to all the ACF values of 10^^, 
in order to correct for a po ssible systematic off set that 
can contaminate the signal (|Blake fc Wal]|[200^ . In this 
case, our constraints are weaker but consistent with the 
previous analysis and we get /nl = 42 ± 30 (Icr CL). 
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Fig. 5. — Dependence on the maximum separation 6 of the error 
bars on the minimal halo mass M^^^ and on /nli calculated with 
the Fisher matrix method (arbitrary normalization). 



We have also checked t hat the correction proposed by 
iWands k, Slosar I ([20091) to account for the infrared di- 
vergence of the non-Gaussian halo correlation function 
is negligibly small for our best-fit /nl value. 

To allow for the "integral constraint" (measurements 
probe the survey mean, and not the ensemble mean), 
we have added to w{9) a constant c and have marginal- 
ized over c, allowing this quantity to vary in the range 
[10~^, lO"**] (the upper limit cannot be larger since this 
is the theoretical variance expected on the scale of the 
survey scales). We find /nl = 58±28 (1 cr), showing that, 
as expected given the large sky fraction covered by the 
NVSS survey, the best-fit value of /nl is only marginally 
affected. 

We use a Fisher matrix approach to estimate which 
scales contribute most to the signal for /nl and 
\ogiQ Mmin from the ACF as a function of 9. Depend- 
ing on where the signal is localized, this may give some 
insights into what systematic effects are important. 

Figure [5] shows the Fisher-predicted error (for the ACF 
of a survey with NVSS characteristics) as a function of 
the maximum angle 6 (the minimum is always set to 
1°). The error-bar normalization is arbitrary and error 
bars of /nl and log^Q M^in at 6* = 8° are set to be the 
same. The error on Mmin stabilizes at ~ 3°, indicating 
that the bias signal is mainly localized at separations 
smaller than 3°. The /nl error decreases rapidly and 
stabilizes at larger 0, indicating that the non- Gaussian 
signal is localized in the ACF at 6* = 2° — 5°. Since the 
mean redshift of NVSS sources brighter than 10 mJy is 
about 1.2, 1° corresponds to a comoving size r ~ 60 Mpc. 
The maximum non-Gaussianity signal thus comes from 
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comoving scales in the range 120 Mpc < r < 300 Mpc, 
while the constraints on Af,„in come primarily from scales 
r< 180 Mpc. 

5. DISCUSSION AND CONCLUSIONS 

While previous analyses exploiting the NVSS to con- 
strain primordial non-Gaussianity of local type have fo- 
cused on the cross-correlation with the WMAP Inter- 
nal Linear Combination map, we have shown that the 
angular correlation function alone is a very sensitive 
non-Gaussianity probe. The key point is that, given 
the redshift distribution of NVSS, which has been re- 
cently observationally determined, and the typical AGN 
halo mass estimated for optically selected QSOs, the 
standard ACDM cosmology with scale-independent bias 
would imply, on scales > 4°, a negative ACF, but, on 
the contrary, it is observed to b e positive. Careful a nal- 
yses of the NVSS sample (e.g. IBlake fc Wall I200I in- 
dicate that systematic offsets that may induce a spuri- 
ous positive signal should be negligible for the sources 
with 5'i.4GHz > lOmJy, selected for the present anal- 
ysis. If so, the NVSS ACF may point at the pres- 
ence of a small primordial non-Gaussianity that adds 
power to the largest scales. After marginalizing over all 
the other parameters, we find 25 < /nl < 117 at the 
95% confidence level, compatible with bounds derived 
by other studies. The minimum halo mass turns out to 
be A/„iin = 10i2i7±o.26;j-i^^^ (la), remarkably close 
to the value found by iCroom et al.l (|2005l) for optically 
selected QSOs. 

We have addressed the significance and the robust- 
ness of our findings by considering different bias mod- 
els and by investigating the impact of gauge effects 
on large scales. Error bars were estimated by a jack- 
knife re-sampling procedure, widely used in the literature 
(|Scranton et al. 1 120031: fXia et al (|2009[) . It is known to 
be robust and accurate for the diagonal elements of the 
covariance matrix but not as widely tested and calibrated 
for the off-diagonal ones, which cannot be neglected be- 
cause n^ighbormg^ACF data points are highly correlated. 
From lScranton et al. I ()2003l ). we infer that jackknife can 
underestimate parameter errors by up to 30%. If the er- 
ror bars were to be increased by this (maximal) amount, 
/nl would become compatible with zero at the ^ 2a con- 
fidence level. We conclude that our work should be seen 
as a "proof of principle" , indicating that future surveys 
probing scales ~ 100 Mpc at substantial redshifts can 
put stringent constraints on primordial non -Gaussianity 
(e.g. ICaripone et an[200l IViel et al. 1 120091) . 
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